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Magnetic field generation in relativistic shocks 
An early end of the exponential Weibel instability in electron-proton plasmas 
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Abstract. We discuss magnetic field generation by the proton Weibel instability in relativistic shocks, a situation that applies to 
the external shocks in the fireball model for Gamma-ray Bursts, and possibly also to internal shocks. Our analytical estimates 
show that the linear phase of the instability ends well before it has converted a significant fraction of the energy in the proton 
beam into magnetic energy: the conversion efficiency is much smaller (of order m c /m p ) in electron-proton plasmas than in pair 
plasmas. We find this estimate by modelling the plasma in the shock transition zone with a waterbag momentum distribution 
for the protons and with a background of hot electrons. 

For ultra-relativistic shocks we find that the wavelength of the most efficient mode for magnetic field generation equals the 
electron skin depth, that the relevant nonlinear stabilization mechanism is magnetic trapping, and that the presence of the hot 
electrons limits the typical magnetic field strength generated by this mode so that it does not depend on the energy content of 
the protons. We conclude that other processes than the linear Weibel instability must convert the free energy of the protons into 
magnetic fields. 
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1. Introduction 

Magnetic field generation in relativistic shocks in a hy- 
drogen (electron-proton) plasma is importan t for the fire - 
ball model for Gamma-ray Bursts jRees & Meszaros 1992). 
This model proposes that the non-thermal radiation we 
observe in the prompt and afterglow emission from the 
Gamma-ray Burst is synchrotron radiation from collision- 
less relativistic shocks. To explain the observed intensity 
of the afterglows as synchrotron emission, the models need 
a magnetic field strength B o f at least ten per cent of 
the equipartition field stre ngth ( Gru zinov & Waxmanlfl999l 
IPanaitescu & KumaJ E"002 ) . This means that the magnetic en- 
ergy density must contribute about one to ten per cent to the 
total energy density of the plasma behind the shock: B 2 /8n ~ 
(0.01 - 1) X e, with e the post-shock thermal energy density. 
In electron-proton plasmas this implies a much stronger mag- 
netic field than in pair plasmas because the energy density in an 
ultra-relativistic shock propagating into a cold medium (where 
nmc 1 » P, with n the number density, m the rest mass and P 
the thermal pressure) is roughly proportional to the rest mass 
m of the particles. In this paper, we look at what happens in a 
plasma consisting of particles of widely different mass. 

In collisionless shocks, plasma instabilities can generate 
magnetic fields. Within the shock transition layer the rela- 
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tive motion of the mixing pre- and post-shock plasma pro- 
duces very anisotropic velocity distributions for all particle 
species concerned. Fluctuating electromagnetic fields deflect 
the incoming charged particles and act as the effective colli- 
sional process need ed to complete the shock transition (e.g., 
iHoshino et alJll992i) . These fluctuating fields occur naturally 
because anisotropic velocity distributions are unstable against 
several plasma instabilities, such as the electrostatic two-stream 
instability and the electromagnetic Weibel instability. The first 
is an instability of (quasi-)longitudinal charge density perturba- 
tions and leads to fluctuating electric fields satisfying \E\ » \B\ 
where E and B are the fluctuating electric and magnetic fields. 
The second is an instability of the advection currents (propor- 
tional to the beam velocity) that result from charge bunching 
in the beams, and leads to spontaneously growing transverse 
waves dWeibell fl959) with \B\ > \E\. In a relativistic shock, 
where the relative velocity of the pre- and post-shock plasma 
approaches the velocity of light, the Weibel instab ility dom- 
inates because it has the largest growth rate dCalifano et alJ 
120021) . 

As both analytical estimates and numerical simula- 
tions show, the Weibel instability in pair plasmas can 
prod uce a magneti c field of near - equip a rtition stre ngth 
dFonseca et all 120031 iHaruki & Sakail 120031 ISazimura et all 
ll998HMedvedev & Loebll999HYang et alJl994 . In numerical 
simulations, the magnetic field generation always undergoes 
an exponentially growing phase that agrees with the estimates 
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from li near analytical th eory, and enters a nonlinear phase af- 
ter that. lYang et al.ldl994l) have shown that in pair plasmas, the 
magnetic field strength reaches its maximum value at the end of 
the linear phase of the instability. The question arises whether 
the same holds true for the Weibel instability operating in an 
electron-proton plasma (Medvedev & Loeb 1999). Numerical 
simulations (Frederikse n et al.ll2004l) show that the nonlinear 
phase may be more important in electron-ion plasmas. 

In this paper we present an analytical estimate that shows 
that the linear phase of the instability ends much earlier for pro- 
ton beams in a hydrogen plasma than for electron(-positron) 
beams in a pair plasma. We do not present a full self-consistent 
shock model: rather we consider the plasma processes that 
could generate a magnetic field in a plasma with properties 
such as one expects near the front of a collisionless relativistic 
shock. The paper is organised as follows. In Sect. |2 we define 
a model for the shock situation in terms of the momentum dis- 
tributions of the particles. We calculate the conditions for the 
instability in Sect. 13. H and the magnetic field strength at the end 
of the linear phase of the instability in Sect. 13.21 In Sect. 13.31 
we compare the energy density associated with this magnetic 
field strength with the total available energy density associated 
with the beams. Section |4] contains the discussion and Sect. |5] 
the conclusions. 

2. A simple model for a relativistic shock transition 

In an electron-proton plasma, the protons dominate the shock 
energetics because they have a much larger rest mass than the 
electrons. Therefore, we will study the proton-driven Weibel 
instability. In this section we present a simple model for the 
plasma in the transition layer at the front of the shock. 

The plasma in an astrophysical relativistic shock does not 
necessarily behave as a single fluid. Coulomb collisions be- 
tween electrons and protons are not sufficiently fast to create 
thermal equilibrium between the protons and electrons. This 
problem of the non-equilibration of the electron and ion ener- 
gies already exists in the much slower (~ 100 km/s) shocks 
associated with Supernova Remnants llDraine&McKeell 19931 
IVinkl2004 . 

We assume that scattering by plasma waves is far more ef- 
ficient for the light electrons than for the heavy ions so that 
when the trajectories of the incoming protons start to become 
significantly perturbed, the electrons have a lready u ndergone 
the fast-growing electron Weibe l instability jFrederiksen et alJ 
l2004 iMedvedev & Loeblfl999i) . which has converted the ki- 
netic energy of their bulk motion into the thermal energy of a 
relativistic ally hot electron plasma with an (almost) isotropic 
thermal velocity distribution. The incoming protons form, seen 
from the rest frame of the hot electrons, a relativistic beam. 
We also assume that part of the protons are reflected further 
downstream although that assumption is not critical for our fi- 
nal conclusions (see Sect.lAl. 

The electron-driven Weibel instability produces a weak 
fluctuating magnetic field with B 2 /8n ~ e e , with e e the energy 
density of the shocked electrons. We ignore this magnetic field 
in the calculations for proton beams, but it could serve as a seed 
perturbation for the proton-driven Weibel instability. 



2.1. The proton velocity distribution 

A simple model for the anisotropic proton velocity distribu- 
tion with in the shock transition layer is a wa terbag distribu- 
tion dSilva et all2002UYoon & Davidsonll987l). W e consi der a 
similar situation as in Fig. 6 of Frederiksen et 

ID £004): we 

take two counter-streaming proton beams moving along the 
x-direction, with a small velocity spread in the z-direction to 
model thermal motions: 

Hp 

F O) = - — [6(p x - pm) + S(p x + p M )\ 

4 />z0 (1) 

X 6(Py) [®(p z + p zQ ) - ®(p z - p z0 )] . 

Here n p is the total proton density, p^> is the bulk momentum of 
the proton beams, p z o is the maximum momentum in thermal 
motions and ®(x) = (1 +x/\x\)/2 is the Heavyside step function. 
The assumption of two beams of equal strength is mathemati- 
cally convenient, but not essential (see Sect.lAl. 

This is a simple model that mimics the properties of non- 
re lativistic coll isionless shocks (see the Micro structure Section 
in Tsurutani & Stone 1985) in which (partial) reflection of the 
ions occurs as a result of deflection by an electrostatic potential 
jump in the shock transition, or by 'overshoots' in the strong 
magnetic field in the wake of the shock. In addition, the wa- 
terbag model accounts for partial ion heating by including a 
velocity dispersion in the direction perpendicular to both the 
beam direction and the wave magnetic field. This direction lies 
along the wave vector of the unstable modes (the z-direction in 
our configuration). 

2.2. The shock conditions for the electrons 

We assume that the electrons have (almost) completed the 
shock transition s o that their properties obey the relativistic 
shock conditions jBlandfordfcMcKeelll97fili . which follow 
from the generally valid conservation laws for particle number, 
energy and momentum. 

Here and below, we will label properties of the post-shock 
electron plasma with subscript 2, and those of the pre-shock 
plasma with subscript 1. We will assume that the pre-shock 
plasma is cold in the sense that e e ,i n e ^m e c 2 . Then the shock 
conditions for the proper density « e and the proper energy den- 
sity e e for the electrons are: 

"e,2 = (4y re i + 3)n e ,i' 

2 ( 2 ) 
<? e ,2 = (4y ie i + 3)y re ira e ,i m eC . 

where y 2 el = 1 + u\. is the Lorentz factor of the relative velocity 
between the pre- and post-shock plasma (with u ^ = p x o/m p c). 

We neglect the dynamical influence of a pre-shock mag- 
netic field on the electron-fluid jump conditions. This influence 
will be small if VA e ,i <Z c where VAe.i = B\/ J4-7m e jin e is the 
Alfven speed based on the electron mass (instead of the pro- 
ton mass). We also neglect any large-scale electrostatic field in 
the shock that might accelerate the electrons to higher energies 
while decelerating the incoming protons. 
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3. The Weibel instability 

3.1. The dispersion relation 

We consider the instability of a purely transverse electromag- 
netic wave with wave vector k - ke z and frequency to. The 
instability will grow with a rate equal to the imaginary part of 
the wave frequency: cr - 3(co). 

The Weibel instability for the model of t he previous sectio n 
obeys a dispersion relation of the form (e.g-. ISilva et all2002l) 



k 2 C 2 -U 2 [l + X xx(0i,k)] = 0, 



(3) 



where Xxx(w, k) = Xxx,e(w, k) + Xxx,r,(w, k) is the ^-component 
of the plasma susceptibility tensor, which contains contribu- 
tions of both the electrons and the protons. 

Since the electrons have a relativistically hot thermal veloc- 
ity distribution their contribution is 



Xxx,e(<°' k ) 



" CO 2 



(4) 



where o> pe is the electron plasma frequency (in Gaussian units): 



For what follows it is convenient to introduce the frequency 
cbpp defined by 



^pp 

Tbo 



(11) 



If the velocity dispersion in the beam is small, which is always 
true for a Weibel-unstable proton distribution (see the end of 
this section), we have yt,o - Trei so that 



^pp 

OJ pe 



4m e 
3ra D 



(12) 



Note that one can get the equations for an electron-po sitron 
beam in an electr on-positron plasma (see also lSilva et alJ2 002: 
lYang et alJl994h by replacing n p and m p with the beam density 
and electron mass respectively. We will use this in what fol- 
lows to compare results for electron-proton plasmas with those 
for electron-positron plasmas. In those cases we assume that 
the density of the electron-positron beams is comparable to the 
density of the background electron-positron plasma. 

Substituting the contributions Q and @ in @ we can 
write the dispersion relation as a biquadratic equation for u>: 



_ 2 4-nq 2 n^ 2 



mJi 



(5) 



with n e the electron proper density, q the electron charge, m t 
the electron mass and h — (e e + P e )/(n e m e c 2 ) 4e e /(3n e m e c 2 ) 
the electron enthalpy per unit rest energy for the relativistically 
hot electrons with e e 3f e » n e m e c 2 . If we assume that the 
electrons are fully shocked, the shock conditions enable us 
to express o> pe in terms of the pre-shock electron number den- 
sity: 

(4y rel + 3 N 



\2nq 2 n^i 



4y re i 



(6) 



The factor between brackets approaches unity in ultra- 
relativistic shocks with y K \ » 1 . 

The proton c ontribution to dispersion relation Q is (see 
ISiivaetall200i 



XxxAu, k) 



TbO w z 



r + 



k 2 v 2 

K V xQ 



■ k 2 v 2 



with 



2v, \c-v z0 



*f0 



(7) 



(8) 



Here w n 



p P - yj47Tq 2 n p /m p is the (non-relativistic) proton 
plasma frequency based on the density in the lab frame, m p 
is the proton rest mass, m,- = pi/(m p c), yt,0 = (1 + u x q + u %) 1 ^ 2 
and v, = Uic/jho- To ensure quasi-neutrality of the plasma we 
must have 

Tip » « e ,2 > (9) 

and the associated plasma frequency is 



with 



where 



- Sco 1 + C = , 



S = k 2 (c 2 + v 2 z0 ) + SI 2 , 



C = k 2 \v 2 zQ {k 2 c 2 + £l 2 )-u 2 pp v 2 J, 



tf = o? pe + cbl p r. 



Since & > the wave is unstable for C < with u> 2 
where the growth rate <x follows from 



cr 2 = 



VS 2 - AC - 



(13) 



(14) 



(15) 



(16) 

-cr 2 <0 



(17) 



For a given set of shock parameters (jbo, v z o, no), the growth 
rate is a function of the wave number (Fig.^. 

Anticipating our results for a proton beam in a back- 
ground of (relativistically) hot electrons, we will assume that 
the growth rate of the unstable modes satisfies 



cr <g kc and cr <K ib. 



pe ■ 



and that the characteristic plasma frequencies satisfy 



Wp P « 



(18) 



(19) 



see Eq. J27I below. 

Under these assumptions we can approximate the solution 
of the dispersion relation with cr 2 - -CIS, which leaves the 
instability criterion (C < 0) unchanged. We will also make the 
approximation of a small beam velocity dispersion: v 2 Q <K c 2 . 
Then the dispersion relation reduces to 



co 2 = (4y re i + 3) 



4nq 2 n 



e,l 



m n 



(10) 



cr 2 = 



kZ (KKu-&vj -k 2 c 2 V 2 ) 
k 2 c 2 + Q 2 



(20) 
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k in units of oj pe cr 1 



1000. 



Fig. 1. Growth rate as a function of wave number for a shock 
with ybo = 1000 and v ; o = 0.001c. The solid line is for proton 
beams in a background of hot electrons; the dashed line is for 
electron-positron beams in an electron-positron plasma. 

For further analysis we introduce the following dimension- 
less quantities: 



kc cr _ <y pp / 

q" ' V ~Q. , a ~ IT \v~q~ 



(21) 



Expressed in these quantities the dispersion relation J20l > reads 



/ v zQ \ 2 K (*™ax K ) 
~\~) 1 +K 2 

\a) lei (i+K 2 )(i + 4 



(22) 







where we have eliminated v z o using the definition of a and we 
define 

*max = V« 2 - 1, (23) 

which is the limiting wavenumber of the Weibel instability: the 
instability condition is satisfied for k < K max . The parameter 
a is a measure for the range of unstable wave numbers. The 
Weibel instability occurs only if a > I. This condition poses a 
restriction on the velocity spread v,n: 



VzO < ^PP 



(24) 



The quantity O depends on v-o through Eqs. (I16> and <|HJ, but a 
quick inspection shows that T and D. are increasing functions 
of v z o so if condition J24i holds, then we also have 



V'zO 
VxO 



'PP 



^PP 



(25) 



It follows that for relativistic shocks with \>xo — * c the proton- 
driven instability requires v z o <sc c so that 



i 



"zO 

3c 2 



« 1, 



and 



D. - LOpe » W pp 



(26) 



(27) 



(see Ea.ll2>. 

The largest growth rate occurs for a mode with wave num- 
ber k„, which follows from (dcr/dk)^^ = 0. Using the approx- 
imated dispersion relation i22t we find that this wavenumber 
equals 

^ = K,(a) = V^H" , (28) 

with the corresponding growth rate 

oik*) , VzO 

—— = v» = (a - 1) — 
12 c 



1 



Q 



v.to 
c 



(29) 



In the last equality we have used the definition of a. For a 
strong, relativistic shock we find o-(k„) w pp . 

In view of this fact and Eq. Q7> . conditions Jl 8I > and dl9l 
automatically hold for the proton-driven Weibel instability. 

3.2. Stabilization of the Weibel instability 

The linear phase of the Weibel instability (during which pertur- 
bations grow exponentially with time) ends when the generated 
electromagnetic fields significantly perturb the trajectories of 
the particles taking part in the instability. Because the magnetic 
fields generated by the instability are inhomogeneous, the parti- 
cles will quiver under the influence of the Lorentz force. When 
the amplitude of these quiver motions exceeds the wavelength 
of the instability, the linear theory breaks down. I Yang et all 
( 1994) give a full treatment of these quiver motions and show 
that this criterion agrees with the magnetic trapping argument, 
which says that the instability will stop when the wave mag- 
netic field becomes so strong that it traps the beam particles. 

The linearized equation describing the quiver motions in 
the z-direction for a beam particle in a magnetic field B(z, t) e y 
reads: 

d% = gv,QB(z,t) 
dt 2 ywnpC 



(30) 



where i; z is the displacement. In the linear stage of the 
instability the wave magnetic field varies as B(z, t) = 
\B\ exp(crf) sin(fe) for a wave with wave number k and growth 
rate <r, so the amplitude of the quiver motion is 



q\B\v x 



(31) 



yborripco-^ 

The trapping criterion k\g z \ < 1 corresponds to \B\ < B tlap with 

yb0'«pccr 2 (£) 



B 



trap 



kvxoq 



(32) 



This corresponds to Eq. (18) o f lYang et alJ lll994l) . Assuming 
that trapping saturates the Weibel instability at all wavelengths 
we get the typical field amplitude as a function of wavenumber 
(Fig.EJ. 

The maximum field amplitude is reached at those wave 
numbers where cr 2 /k has the maximum value. For dispersion 
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saturation magnetic field 



as 



0.001 




0.01 



0.001 



0.1 1. 
k in units of < 



1000. 



Fig. 2. Saturation magnetic field as a function of wave num- 
ber for a proton-driven instability and for an instability in 
an electron-positron plasma, using the same parameters as in 
Fig.^with n e \ - 1 cm -3 . 

relation d22t this maximum is reached at a wavenumber A:' that 
follows from 



= k (a) 



- VsTo 2 

2 



1 " T 



(33) 



If the instability is strong, so that a » 1, we find that this 
wavenumber and the associated growth rate are 



k c „ j 



V2 c 



LD 



V2Q 



^. (34) 
c 



Note that the trapping criterion predicts the largest field ampli- 
tude at a wavenumber fc 1 , which is not equal to the wavenum- 
ber k m of the fastest growing mode (Eg. I28l i: k* - Q./c «: 
K — -\/aQ./c for a » 1. Eventually, this slower growing 
long-wavelength mode reaches a higher magnetic field than 
the fa ster growing mode (see also the discussion in lYan^ et alJ 
11994 . If a » 1 the saturation field strength due to trapping at 
a Q/c is 



B 



trap 



PP 

2Q. 



q 

2yu)m e LbpeC 



(35) 



where we have used Eqs. (I12> and (I27l i to eliminate u> pp and Q. 

For parameters typical for the external shock associated 
with Gamma-ray Bursts we have 



Strap - 3.7 VnTyiooo Gauss, 



(36) 



with n\ = «o/(l cm 3 ) and y 1 nnn - Tho / 1 000- 

We should note that Medvedev & Loeb ll 19991) use a differ- 
ent method for estimating B: they propose that the instability 
saturates when the beam ions become magnetized. This hap- 
pens when the Larmor radius rt = yvmc 2 IqB of the beam par- 
ticles in the generated magnetic field becomes smaller than the 
wavelength of the fastest growing mode of the instability. This 
criterion k t r^ < 1 corresponds to |B| < B magn with 



ybov.xom ? ck t 



(37) 



with K given by Eq. (I28> . However, the trapping argument pre- 
dicts the smallest saturation amplitude \B\. In particular, we 
have 

(38) 



Strap cr 2 (£t) 



B 



magrj 



where 



O(a) = 



iff ( a ) 



v() 



= 1 <■><«>• 



a 2 Va - 1 



3a- V8 + a 2 



V8 



+ a z 



(39) 



To derive this relation we have used definition J2 1 i to write 
Vjo/vxO = w pp /(aQ). For a proton beam in a hot electron back- 
ground we have w 2 p /0 2 <§c 1 and <t>(a) < 0.3 for all a > 1. For 
a strong instability with a » 1 we have <I>(a) (4a) -1 ' 2 <K 1. 
Therefore, trapping occurs well before the field can totally 
magnetize a proton beam with a density comparable to the den- 
sity of the hot background electrons. In view of this we will use 
B trdp as an estimate of the saturation magnetic field strength. 

The criteria fl32i and J37l > predict the typical amplitude of 
the magnetic field as one particular wave mode k saturates. In 
a realistic situation the instability will involve a superposition 
of wave modes and one should interpret B as the amplitude 
that follows from the power spectrum Is(k) of the field fluc- 
tuations: B 2 l%n ~ kls(k) with k » k* for magnetization and 
k m fef for trapping. The total magnetic energy in the unstable 
modes is 

b 2 r*- 

Ub = — = dkI B (k). (40) 



3.3. The equipartition parameter 

A measure of the strength of the magnetic field is the equiparti- 
tion parameter, which compares the energy density in the mag- 
netic field with the total energy density. 

The protons dominate the energy budget, and the total 
available energy density is 

e p = j dpF(p)y(p)m p c 2 - y h0 n p m p c 2 , (41) 

where the approximation is valid for v z o <K v^. 
We define the proton equipartition parameter as 



B 



trap 



(42) 



Using Eq. J35i for the magnetic field, with the definition dl It 
for o) pp and the approximation v A -o =s c for relativistic shocks 
we get 

e B = ■ (43) 



Then from Eqs. (I12t and Mil we have 



6m D 



10" 



(44) 
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4. Discussion 

The small value of the equipartition parameter J44l> implies that 
the proton-driven Weibel instability in a background of rel- 
ativistically hot electrons saturates long before the magnetic 
field reaches equipartition with the available proton free en- 
ergy. Equation d25l > demonstrates that the proton-driven Weibel 
instability in a hydrogen plasma is not a suitable candidate for 
the mechanism responsible for 'thermalization' of the incom- 
ing protons in the shock layer. In electron(-positron) beams in 
a hot pair background the instability condition, a > 1, allows 
a large beam velocity dispersion: v z q/v x o < 1 if the number 
densities in the beam an d in the hot b ackground are of sim- 
ilar magnitude (see also lYoon & DavidsorJll987l) . Therefore, 
the Weibel instability is in principle capable of randomizing 
a significant fraction of the beam momentum of an electron(- 
positron) beam, as asserted in Sect. |5J whereas this is not true 
for a proton beam in a hydrogen plasma. 

In the limit of a cold, relativistic proton beam the properties 
of the electrons and not those of the protons determine many 
of the results. The electron plasma frequency determines the 
wave number fl33i of the mode with the maximum magnetic 
field so that the electron skin depth c/w pe sets the length scale 
of the dominant mode. This happens because the low inertia of 
the electrons makes them very responsive to the perturbations 
of the protons. The dispersion relation for the Weibel modes 
(Fig-E> also supports this view: the plateau around the maxi- 
mum growth rate starts roughly at a wave number k ~ u> pss /c. 
Studies that do not include the response of the background elec- 
trons (by treating the protons as an isolated system) miss this 
point. The peak magnetic field (T35> does not contain any pa- 
rameters connected with the protons. Therefore, proton beams 
in a hydrogen plasma generate nearly the same magnetic field 
strength as electron(-positron) beams in an electron-positron 
plasma (Fig. 13 despite the larger kinetic energy of the protons. 

iGruzinovl 11200 11) anticipated this when he excluded the case 
where a small parameter in the theory might be important in his 
analysis of the Weibel instability: our analysis shows that the 
relevant small parameter is a>^ p /Q. 2 ~ 4m e /3m p . The result is a 
small equipartition parameter J44i . I n this respect our result is 
similar to the one found by Sagdeev ( 1966, p. 88), who argued 
for the Weibel instability in a non-relativistic plasma that the 
electrons have a quenching effect on the ion instability. 

In our analysis we have excluded electrostatic waves, 
which cou ld also play an importan t role in the shock transi- 
tion zone JSchlickeiser et alJl2002i) . In that case electrostatic 
Bremsstrahlung could be an alternative explanation for the 
Gamma-ray Burst afterglow emission lISchlickeiserl EoO^. re- 
laxing the need in synchrotron models for a high magnetic field 
strength. 

5. Conclusions 

We have presented an analytical estimate of the magnetic field 
produced at the end of the linear phase of the Weibel instabil- 
ity at the front of an ultra-relativistic shock propagating into a 
cold hydrogen plasma, a situation that applies to the external 
shocks that produce gamma-ray burst afterglows in the fireball 



model jRees & Meszaroslll992l) . The magnetic field strength 
that we find is too weak to explain th e observed synchrotron ra- 
diation JGruzinov & Waxmanl 19991) : the equipartition parame- 
ter (~Eq. l44t is at least two orders of magnitude too small. This 
is radically different from the results for the Weibel instability 
in pair plasmas. The reason is that the contribution of the elec- 
trons to the electromagnetic response of the plasma inhibits the 
instability of the protons. 

The saturation magnetic field d35l > which this low equiparti- 
tion parameter corresponds to is the magnetic field at the point 
where the linear approximation breaks down and where non- 
linear trapping effects start to limit the growth of the unstable 
Weibel mode. After this happens, it is likely that the instabil- 
ity enters a nonlinear phase or that anothe r type of i nstability 
takes over: numerical simulations ( Frederiksen et al. 2004J) of 
similar plasmas show near-equipartition magnetic fields behind 
the Weibel-unstable region in the shock transition. The nonlin- 
ear phase would then be the dominant phase in electron-ion 
plasmas and deserves further study to determine the physical 
mechanism and the properties of the resulting magnetic field. 
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Appendix A: Asymmetric beams 

We consider asymmetric proton beams and show that we can 
neglect the effect of the asymmetry for typical parameters. We 
replace the proton distribution Q by 



n p 1 + A 1 - A 

F(j>) = - — 8(p x - pxo) + — — 6(p x 

2p z [ 2 2 

X S(py) [®(p z + p z0 ) - ®(p z - p z0 )] , 



Pxo) 



(A.l) 



with < A < 1 the parameter measuring the asymmetry be- 
tween the two beams. This a symmetry changes th e dispersion 
relation of the waves to (e.g., Akhiez er et al.1 19751) : 



1 +Xxx(oj, k) 



1 +Xu(o), k) 



0. (A.2) 



proton— driven weibel waterbag disp. relation 

A = 0, 0.5, 0.9 and 1 
C e = 0/V3 
7 rel = 1000 - ^o= 001 




0.01 0.1 1 10 100 

wavenumber tc = kc/n 

Fig. A.l. The growth rate as a function of wavenumber for 
asymmetric proton beams in a relativistic ally hot electron back- 
ground with sound speed C e — c/ y3. Different lines corre- 
spond to different values of the parameter A, which measures 
the asymmetry between the beams. 



The extra term oc x x ~ appears in relation (IA.2> because the 
waves are no longer transverse: the charge bunching of the 
asymmetric beams produces a net charge density, leading to 
a component of the wave electric field along the wave vector. 

We will give the extra components of the susceptibility ten- 
sor and publish a derivation elsewhere (Achterberg & Wiersma, 
in preparation). The proton contribution to Xxx(oJ, k) is the 
same as in the symmetric case. For a thermal electron back- 
ground with an isotropic momentum distribution, only the two 
proton beams contribute to the off-diagonal components of the 
susceptibility tensor so that 



Xxz(u, k) = Xxz, P (u, k) = 



A -pp 



kv xQ 



u> 2 - k 2 v 2 Q o) 



(A3) 



This off-diagonal component vanishes in the symmetric case 
(A = 0). The longitudinal response of the beam-plasma system 
is contained in 



1 +Xz.z(w, k) = 1 



^pe 



"pP 



CO 2 - k 2 Cl 0? - k 2 v 2 zQ 



(A.4) 



Here C 2 = 4P e /3n e m e h is the effective sound speed in the 
hot electron gas. In the ultra-relativistic case one has C e - 
c/ y3. For proton beams in a hot electron background one has 
cb 2 e » &p P and \a)\ < oj pp <sc kC e for the Weibel instability (see 
Sect. 13. 1> . This implies that 



1 +Xzz - 1 + 



pp 



k 2 Cl 



o 2 — k 2 v 



.2 • 

.-0 



(A.5) 



One finds 

(Z-Zi)(Z + Zi)-A 2 ZxZ2 = 0. 

Here Z 1,2 ar e defined by 

k 2 C 2 



k 2 Cl + 



2 , r,2 



Zi 



k 2 v 2 

K V x0 

k 2 c 2 + O 2 



(A.7) 



(A.8) 



with Zi > Z\ for the parameters considered here. The solution 
of this quadratic equation, 



_ Z\-Zi (Zi + Zz\ -r 
Z± = + a - A Z\Zi ■ 



determines the frequency through 



- 2 -r 
w pp <-± 



,2 2 
k V, n 



(A.9) 



(A. 10) 



The Weibel-unstable branch corresponds to the solution branch 
w_ as Z- < 0. In the symmetric case (A = 0) one has 
Z- - -Z2 which follows from iA.l\ . and one recovers disper- 
sion relation i20\ . The stable branch u> + is a modified (largely 
electrostatic) ion-acoustic wave. 

Although the asymmetry decreases the range of unstable 
wave numbers (Fig. IA. 1 1 and lowers the growth rate with re- 
spect to the symmetric case A = 0, the change is small unless 
A a; 1: the case where there is almost no reflection. For a single 
beam, A = 1, one has Z- — Z\ - Zi, which gives the follow- 
ing dispersion relation for v 2 = -a> 2 /O 2 in the ultra-relativistic 
limit with C e — cj V3 and D. 2 w w 2 e : 



The resulting dispersion relation can be written in terms of 
the dimensionless variable 



Z(6>, k) = 



oj 2 - £ 2 v 2 



'pp 



(A.6) 



""pp 
D 2 



Civ 

c 



1 



K 2 



V 2 ^ 

1 _ 

c 2 



(3 + K 2 )(l + K 2 ) 



-K 2 ^. (A.ll) 



